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Abstract 

The loop-induced processes gg — > h,H,A provide the dominant Higgs boson pro- 
duction mechanisms at the Tevatron and LHC in a large range of the minimal 
supersymmetric extension of the Standard Model. For squark masses below ~ 400 
GeV squark loop contributions become important in addition to the top and bottom 
quark loops. The next-to-leading order QCD corrections to the squark contribu- 
tions of these processes are determined including the full squark and Higgs mass 
dependences. They turn out to be large and thus important for the Tevatron and 
LHC experiments. Squark mass effects of the K factors can be of 0(20%). In addi- 
tion we derive the QCD corrections to the squark contributions of the rare photonic 
Higgs decays h,H^ 77, which play a role for the Higgs searches at the LHC. 



* Supported in part by the Swiss Bundesamt fur Bildung und Wissenschaft. 



1 Introduction 



Higgs boson pQ searches belong to the major endeavors at present and future colliders 
within the Standard Model (SM) and its minimal supersymmetric extension (MSSM). 
In the MSSM two isospin Higgs doublets are introduced in order to generate masses of 
up- and down-type fermions [2J. After electroweak symmetry breaking three of the eight 
degrees of freedom are absorbed by the Z and W gauge bosons, leaving five states as 
elementary Higgs particles. These consist of two CP-even neutral (scalar) particles h, H, 
one CP-odd neutral (pseudoscalar) particle A and two charged bosons H ± . At leading 
order the MSSM Higgs sector is fixed by two independent input parameters which are 
usually chosen as the pseudoscalar Higgs mass M4 and tg/3 = v 2 /vi, the ratio of the two 
vacuum expectation values. Including the one-loop and dominant two-loop corrections 
the upper bound of the light scalar Higgs mass is M h < 140 GeV [3]. The couplings of 
the various neutral Higgs bosons to fermions and gauge bosons, normalized to the SM 
Higgs couplings, are listed in Table HJ where the angle a denotes the mixing angle of the 
scalar Higgs bosons h, H. An important property of the bottom Yukawa couplings is their 
enhancement for large values of tg/3, while the top Yukawa couplings are suppressed for 
large tg/3 pE]. 
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Table 1: Higgs couplings in the MSSM to fermions and gauge bosons [V = W, Z] relative 
to the SM couplings. 



For this work we need the Higgs couplings to stop and sbottom squarks in addition. 
The scalar superpartners ji jj of the left- and right-handed fermion components mix with 
each other. The mass eigenstates fi t 2 of the sfermions are related to the current eigenstates 
f L>R by mixing angles 8 f , 



fi = IlcosOj + f R sm9 f 
h = -/isinfl/ + f R cos9 f . 

The mass matrix of the sfermions in the left-right-basis is given by p)Ji 



M: 



Ml 



+ m| rrif(Af — fj,rf) 



m 



Mj +m 2 f 

JR. J 



For simplicity, the D-terms have been absorbed in the sfermion-mass parameters M 2 } 

JL/R 



(2) 



1 



with the parameters = l/r u = tg/5 for down- and up- type fermions. The parameters 
Af denote the trilinear scalar couplings of the soft supersymmetry breaking part of the 
Lagrangian, \x the Higgsino mass parameter and rrif the fermion mass. The mixing angles 
acquire the form 
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and the masses of the squark mass eigenstates are given by 
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Since the mixing angles are proportional to the masses of the ordinary fermions, mixing 
effects are only important for the third-generation sfermions. The neutral Higgs couplings 
to sfermions read p] 
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with the couplings gf listed in Table [2j I 3 j denotes the third component of the electroweak 
isospin, e/ the electric charge of the fermion /, 9 W the Weinberg angle and M z the Z- 
boson mass. All these couplings have to be rotated to the mass eigenstates by the mixing 
angle Of. 



f 





fx 


<f> 

92 


9s 


94 




h 


cos a/ sin (3 


— sin(a + (3) 


— sin a/ sin (3 


cos a/ sin f3 


u 


H 


sin a 1 sin f3 


cos(a + (3) 


cos a/ sin (3 


sin a 1 sin /3 




A 








-1 


l/tg/3 




h 


— sin a/ cos (3 


— sin(a + 0) 


cos a 1 cos (3 


— sin a/ cos /3 


d 


H 


cos a / cos (3 


cos(a + (3) 


sin a / cos (3 


cos a/ cos /3 




A 








-1 


tg/3 



Table 2: Coefficients of the neutral MSSM Higgs couplings to up- and down-type sfermion 
pairs. 



At hadron colliders as the Tevatron and LHC neutral Higgs bosons are copiously 
produced by the gluon fusion processes gg — > h/H/A, which are mediated by top and 
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bottom quark loops as well as stop and sbottom loops in the MSSM (see Fig. [T]). Due 
to the large size of the top Yukawa couplings gluon fusion comprises the dominant Higgs 
boson production mechanism for small and moderate tg/3. For large tg/3 the leading role 
is taken over by Higgs radiation off bottom quarks due to the strongly enhanced bottom 
Yukawa couplings [7]. 




Figure 1: MSSM Higgs boson production via gluon fusion mediated by top- and bottom 
quark as well as stop and sbottom loops at leading order. 



The QCD corrections to the top and bottom quark loops are known since a long time 
including the full Higgs and quark mass dependences [8]. They increase the cross sections 
by up to about 100%. The limit of very heavy top quarks provides an approximation 
within ~ 20 - 30% for tg/3 < 5 [9]. In this limit the next-to-leading order (NLO) QCD 
corrections have been calculated before [TU] and later the next-to-next-to-leading order 
(NNLO) QCD corrections [TT]. The NNLO corrections lead to a further moderate increase 
of the cross section by ~ 20 — 30%, so that the dominant part are the NLO contributions. 
Very recently an estimate of the next-to-next-to-next-to-leading order (NNNLO) effects 
has been obtained [12J indicating improved perturbative convergence. A full massive 
NNLO calculation is not available so far, so that the NNLO results can only be used for 
small and moderate tg/3. On the other hand the NLO corrections to the squark loops are 
only known in the limit of heavy squarks [T3|. Moreover, the full SUSY-QCD corrections 
have been obtained recently for heavy SUSY particle masses [H]. NLO computations for 
the last two contributions including the full mass dependences are missing so far. This 
work presents the pure QCD corrections to the squark loops including the full squark and 
Higgs mass dependences as a first step towards a full NLO SUSY-QCD calculation. 

The reverse processes, gluonic Higgs decays, play a role for the decay profiles of the 
Higgs particles. Their partial widths can be measured at a future linear e + e~ collider with 
energy up to about 1 TeV [15]. The gluonic branching ratios can reach a level of O(10%) 
in some regions of the MSSM parameter space [H[71TT6] . Analogously to the gluon fusion 
processes the QCD corrections to the top and bottom quark loops have been calculated 
including the full quark and Higgs mass dependences [8]. The NLO corrections in the 
limit of very heavy top quarks can be found in [TOUT]. This limit provides a reasonable 
approximation for small and moderate values of tg/3. The NNLO QCD corrections are 
known for very heavy top quarks, too [18]. Very recently, even the NNNLO QCD correc- 
tions have been obtained in the heavy top mass limit [TS]. While the NLO corrections 
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enhance the partial decay width by up to about 70%, the NNLO corrections turn out 
to be of moderate size 0(20%) supplemented by the NNNLO corrections in the per cent 
range. The quark mass effects at NNLO and beyond are unknown. The SUSY-QCD 
corrections to the gluonic decay mode are identical to the gluon fusion process. The full 
mass dependences of the NLO SUSY-QCD corrections are unknown. This work provides 
a first step towards this aim by determining the pure QCD corrections to the squark loops. 

A subsample of the diagrams describing the Higgs couplings to gluons determines the 
Higgs couplings to a pair of photons. Since at NLO there is no gluon radiation, the LO and 
NLO results for the photonic Higgs couplings are valid for the Higgs decays into photons, 
which are important Higgs decay modes for the Higgs search at the LHC Moreover, they 
determine Higgs boson production at a future linear photon collider, a satellite mode of 
a linear e + e~ collider built up by Compton back-scattered laser light off the e~ beams 
[2"0~f2"Tj . The NLO QCD corrections to the quark loops are known in the limit of a very 
heavy top quark [22|23] as well as in the fully massive case [22], while the NNLO [21] and 
NLO SUSY-QCD [25] corrections are only known for very heavy top and SUSY-particle 
masses. The results presented in this paper comprise a first step towards a full massive 
SUSY-QCD calculation at NLO by means of the QCD corrections to the squark loops. 

This paper is organized as follows. In Section 2 we describe the NLO QCD calculation 
of the squark loop contributions to the photonic and gluonic MSSM Higgs couplings and 
present the results. Section 3 summarizes and concludes. 



2 NLO QCD Corrections 



In order to compute the pure QCD corrections to the squark loops we need a modifica- 
tion of the MSSM interaction Lagrangian allowing to separate gluon and gluino exchange 
contributions in a renormalizeable way. This is achieved by starting from the basic La- 
grangian^] 
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with the covariant derivative _D M = <9 M + ig s G a u T a + ieA^Q. Here G* denotes the gluon 



field strength tensor and G° the gluon field accompanied by the color SU (3) generators 
T a (a = 1, . . . , 8), while F^ u is the photon field strength tensor and the photon field 
associated by the electric charge operator Q. The Higgs field H represents generically 
either the light scalar h or the heavy scalar H Higgs boson of the MSShfl In this work 
the couplings Qq and are not renormalized, thus leading to a renormalizeable model 



2 Quartic selfinteractions among the squarks are not taken into account in this work. 
3 Since there are no squark loop contributions to the pseudoscalar Higgs boson couplings to photons 
and gluons at leading order (LO), in this paper we will only deal with the scalar Higgs bosons h, H. 
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with strongly interacting scalars Q. Only if gluino exchange contributions were taken into 
account, the full MSSM renormalization would have to be performed. It should be noted 
that the gluino contributions can be separated from the QCD corrections induced by light 
particles due to the gluino mass. Gluino corrections are expected to be small [H]. 

The numerical results in this work will be presented for the gluophobic Higgs scenario 
[26] . which develops strong interference effects between quarks and squarks. It is defined 
by the following choices of the MSSM parameters [m t = 174.3 GeV], 

M SUSY = 350 GeV 

fj, = M 2 = 300 GeV 

X t = -770 GeV 

A b = A t 

m-g = 500 GeV , (7) 

where X t = A t — ///tg/3. We have used the program HDECAY [16] for the numerical 
determination of the SUSY particle masses and couplings. In this scenario the squark 
masses amount to 

m- ti = 156 GeV 
m h = 517 GeV 
mr = 346 GeV 
m- b2 = 358 GeV 

m h = 155 GeV 

m h = 516 GeV 

m~ bi = 314 GeV 

m- b2 = 388 GeV . (8) 

The gluophobic scenario maximizes the destructive interference effects between top and 
stop loop contributions to the light scalar Higgs coupling to gluons. The results of this 
work look similar, whenever the squark masses turn out to be of the order of the top mass, 
or the Higgs mass reaches values beyond the corresponding squark- ant isquark threshold. 



tg/3 = 30 : 



2.1 Scalar Higgs Boson Couplings to Photons 



At leading order (LO) the photonic MSSM Higgs couplings are mediated by top and 
bottom quark loops as well as W loops for the scalar Higgs particles h, H. If the stop and 
sbottom masses range below ~ 400 GeV, there are significant contributions from squark 
loops, too (see Fig. [2]). The LO photonic decay widths are given by [7f8|27] 
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Figure 2: MSSM Higgs boson couplings to photon pairs mediated by top- and bottom quark, 
stop and sbottom as well as W boson loops at leading order. 
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where we neglected contributions from charginos, charged Higgs bosons and the charged 
sleptons. The full expressions can be found e.g. in [7] and have been used in our work. 
N c f(N c f) denote the color factors and e/(ej) the electric charges of the (s)fermions in 
units of the proton charge, while the scaling variables are defined as = 4mf /M^ H . For 
large loop particle masses the form factors approach constant values, 
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The reverse processes, 77 — > h/H, play an important role for the MSSM Higgs search 
at a photon collider at high energies [21~ll28] . Such a photon collider can be realized by 
Compton back-scattering of laser light from high-energy electron beams in a linear e + e~ 
collider [20]. In this collider 77 cm. energies up to about 80% of the corresponding e + e~ 
cm. energy can be reached. The s-channel Higgs boson production cross sections are then 
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given by 
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where dC 11 / dr h / H denotes the differential 77 luminosity for r h / H = M%, H /s 7J with s. 
being the squared cm. energy of the photon collider. This relation between the photon 
fusion cross section and the photonic Higgs decay width holds in NLO QCD, too, since 
single gluon radiation vanishes due to the conservation of color charges as well as due to 
the Furry theorem. 
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Figure 3: Generic diagrams for the NLO QCD corrections to the squark contributions to 
the photonic Higgs couplings. 



The NLO QCD corrections to the photonic Higgs decay modes and the photon fusion 
cross section require the calculation of the two-loop diagrams depicted in Fig. [31 We have 
reduced the 5-dimensional Feynman parameter integrals to one-dimensional ones which 
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have been integrated numerically. In a second calculation the QCD corrections have 
been obtained purely numerically. Both results agree within integration errors. The QCD 
corrections can be parametrized by shifts of the quark and squark amplitudes according 
to 



AT(r Q ) - A"I»{ TQ ) 



:n) 



The QCD corrections to the quark loops can be found in Refs. [22J. The corresponding 

h I FT 

coefficient Cq (tq) of the QCD corrections to the squark loops for finite Higgs and squark 
masses is presented in Fig. H] as a function of tq. In order to improve the perturbative 
behavior of the squark loop contributions they should be expressed preferably in terms of 
the running squark masses 771,5 (M^/2), which are related to the pole masses Mq via 




m (/i) = Mq ( — ^ ] (12) 

where (3q = 33 — 2Np with Np = 5 light flavors. Their scale is identified with \i = M^/2 
within the photonic decay mode. These definitions imply a proper definition of the QQ 
thresholds M-h = 2Mq, without artificial displacements due to finite shifts between the 
pole and running squark masses, as is the case for the running MS masses. We have taken 
into account the LO scale dependence of the squark masses generated by light particle 
contributions. 

The coefficient Cq H (tq) is real below the QQ threshold and complex above. For 
Tq — 1 it diverges, so that within a margin of a few GeV around the threshold the 
perturbative analysis cannot be applied. This singular behavior originates from a Coulomb 
singularity at the threshold, since QQ pairs can form ++ states. This implies that 
the imaginary part 3m Cq H develops a step due to Coulombic gluon exchange, thus 
resulting in a singular behavior of the real part at threshold. The singular behavior can 
be quantified. The lowest order form factor of the squark loops is given by 

a Wh,lo {t) = _ r[1 _ r/(r)] ^ ^ _ 1 + in(] ioTT ^ 1 (13) 



where (3 = a/1 — T denotes the squark velocity above threshold. The QCD corrections 
to the imaginary part at threshold can be derived from the Sommerfeld rescattering 
correction [30] . 

Q^ = T ^ P F = l + § + (« forZ = ^. (14) 



'Analytical results can be derived, too, as in the case of the quark loops [29] . 
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Figure 4: Real and imaginary parts of the QCD correction factor to the squark amplitudes 
of the two-photon couplings for the scalar MSSM Higgs bosons. The renormalization scale 
of the running squark mass is taken to be \x = 



A h J H -> A^ H ' LO + _iog( r r 1 -l)+i7r + const , (16) 



The QCD corrected imaginary part of the h/Hjj couplings at threshold is now given as 

27T 2 

3m A^ H = it/3C cou i = irf3+ ya s , (15) 

so that the real part can be derived from a once-subtracted dispersion relation yielding 
the following expression in the threshold region 

27ra s 

where the smooth non-singular constant is not relevant for the divergent behavior. Finally 
the singular behavior of the QCD correction factor at threshold can be cast into the form 

SeC « /H -3^4) l0 8( T « 1 - 1 ) + 0OTM * 

* mC r ^ 3^4) Kl4 ' 09 ' (17) 

The total size of the imaginary part and the singular behavior of the real part agree with 
the numerical results depicted in Fig. HI However, the singular behaviour at threshold is 
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unphysical. A proper inclusion of the finite decay widths of the virtual squarks as well 
as a resummation of the Coulomb singularities [21] will regularize the divergences and 
improve the perturbative results at threshold. 

In the limit of heavy squark masses the coefficient C^ H (tq) approaches a finite and 
constant valued 



h/H, 
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This asymptotic value can also be derived from low-energy theorems 
are based on the relation between the matrix elements with and without a light external 
Higgs boson. This allows the derivation of an effective Lagrangian for the Higgs couplings 
to photons mediated by squark loops in the heavy squark mass limffl 
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where /3n(a)/a = (a/2n)[l + 4a s /ir + ■ ■ ■] denotes the heavy squark Q contribution to the 
QED f3 function and j m . (a s ) = a s /ir + • • • the anomalous squark mass dimension. The 
NLO expansion of the effective Lagrangian reads as 
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which agrees with the - value of Eq. ( [18]) in the heavy squark limit. 

The numerical results are presented in Figs. [5H7] for the gluophobic Higgs scenario^. 
Fig. [5] shows the partial photonic widths at NLO of the scalar Higgs bosons h, H for 
two values of tg/? = 3, 30. The dotted lines exhibit the photonic widths without the 
contributions of SUSY particles, the dashed lines include all SUSY particle contributions 
except those of the squarks, while the full curves correspond to the full partial decay 
widths. For small as well as for large values of tg/5 the photonic widths range between 
about 0.1 and 10 keV. The importance of the SUSY particles and in particular the squark 
contributions is clearly visible from the comparison of the three different curves. The 
kinks, bumps and spikes correspond to the WW, tjti, ti, bibi, 7yfi, f 2 f 2 and 6 2 &2 thresholds 
in consecutive ordeij^. Due to the Coulomb singularities the stop and sbottom thresholds 
develop spikes corresponding to the logarithmic singularities according to Eq. (TTT1) . 

The relative QCD corrections to the photonic Higgs decay widths are presented in 
Fig-Elfor the two cases, in which SUSY particles have been taken into account or not. The 
QCD corrections reach a size of 10-20% for moderate and large Higgs masses apart from 



5 This value differs from the heavy squark limit obtained in Ref. [25]. The difference can be traced 
back to a wrong expression for the anomalous squark mass dimension used in [25j . 

6 The anomalous dimension of the kinetic photon operator F^ V F^ V does not contribute at NLO |9|17j . 

7 The massive QCD corrections to the quark and squark loops have been implemented in the program 
HDECAY PS]. 

8 In the gluophobic scenario the f masses amount to = 348 GeV, = 356 GeV for tg/3 = 3 and 
m fl = 327 GeV, m f2 = 377 GeV for tg/3 = 30. 
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Figure 5: QCD corrected partial decay widths of the scalar MSSM Higgs bosons to two 
photons as functions of the corresponding Higgs masses for tg/3 = 3 and 30. The full 
curves include all loop contributions, in the dashed lines the squark contributions are 
omitted and in the dotted lines all SUSY particle loops are neglected. The kinks, bumps 
and spikes correspond to the WW, tit\,tt, b\bi, f\f\, f 2 f 2 and 62^2 thresholds in consecutive 
order with rising Higgs mass. The renormalization scale of the running quark and squark 
masses is chosen to be Mh/n/2, while the scale of a s is taken to be the corresponding 
Higgs mass. 
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Figure 6: Relative QCD corrections to the partial decay widths of the scalar MSSM Higgs 
bosons to two photons as functions of the corresponding Higgs masses for tg/3 = 3 and 30. 
The full curves include all loop contributions while in the dashed lines thesquark contribu- 
tions are omitted. The kinks and spikes correspond to the WW, t±ti, ti, b\bi, fif i , f 2 r 2 and 
b 2 b 2 thresholds in consecutive order with rising Higgs mass. The renormalization scale 
of the running quark and squark masses is chosen to be Mh/u/2, while the scale of a s is 
taken to be the corresponding Higgs mass. 
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Figure 7: Ratio of the QCD corrected partial decay widths of the scalar MSSM Higgs bosons 
to two photons including the full squark mass dependence and those obtained by taking 
the relative QCD corrections to the squark loops in the heavy mass limit as functions of 
the corresponding Higgs masses for tg/3 = 3 and 30. The kinks and spikes correspond to 
the WW,titi,ti,bibi,Tif 1 ,f 2 T 2 and b 2 b 2 thresholds in consecutive order with rising Higgs 
mass. The renormalization scale of the running quark and squark masses is chosen to be 
Mh/ii/2, while the scale of a s is taken to be the corresponding Higgs mass. 
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the threshold regions, where the perturbative results are unrealiable due to the Coulomb 
singularities. Since at a 77 collider the photon fusion cross section can be measured with 
an accuracy of a few per cent, these corrections have to be taken into account properly. 
The size of the QCD corrections with and without SUSY particle loops is of the same 
order of magnitude, but they can be of opposite sign. 

In order to quantify the size of squark mass effects beyond the heavy squark mass 
limit in the relative QCD corrections, the ratio between the fully massive photonic decay 
width Y{h/H — > 77) at NLO of Eqs. (0 EE]) and the approximate width T^, where the 
heavy squark limit of Eq. (lisp has been inserted for the coefficient C^ H instead of the 
fully massive result, is shown in Fig. [7J (The full squark mass dependence of the LO 
form factors has been taken into account in both expressions.) The squark mass effects 
reach a size of up to ~ 30%, thus underlining the relevance of including the full mass 
dependences. 



2.2 Gluonic Scalar Higgs Decays 



The gluonic Higgs boson decays h/H — ► gg are mediated by quark and squark triangle 
loops (see Fig. [p. The lowest order decay widths of the scalar MSSM Higgs boson decays 
are given by 
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using the same notation as for the photonic Higgs boson decays. For large loop particle 
masses the form factors approach constant values, 
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The squark contributions become significant for squark masses below about 400 GeV. 

The NLO QCD corrections to the gluonic decay widths are known to be large [SH7J, 
[TSITTU] . They can be decomposed into the two- loop virtual corrections and the one- loop 
real corrections determined by the real radiation processes 

h/H->ggg and gqq. 

The generic Feynman diagrams for the squark loop contributions are depicted in Figs. [S] 
and M As in the photonic case the Feynman integrals of the virtual corrections have been 
reduced to one-dimensional integrals, which have been evaluated numerically. In a second 
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Figure 8: Generic diagrams for the virtual NLO QCD corrections to the squark contribu- 
tions to the gluonic Higgs couplings. 



15 



^oinr- g 

s I 

h/H < x t,6 ro-Q-c^ g h/H < 

s ' 

x HJJJLQ^ g 



h/H- 




-< 



t,b 
-^slqjls g 

K 

t,b 
s Ulqjl> g 



g 



h/H- 



h/H- 





Figure 9: Typical diagrams for the real NLO QCD corrections to the squark contributions 
to the gluonic Higgs decays. 



calculation the results have been obtained purely numerically. Both results are in mutual 
agreement. 

The strong coupling constant a s is renormalized in the MS scheme, with the top quark 
and squark contributions decoupled from the scale dependence, and the quark and squark 
masses are renormalized on-shell. The result can be cast into the form 



T{h/H^gg) = T LO {h/H - gg) jl + E h ' H ^ 

I 7T 



(22) 



E h / H 



T-n h/H A h/H (r 



+ 



4 6^2 \E Q g h Q /H < /H (r Q ) + E,g% H A h f(T,) 

33 - 2N F j4 AEh/H 



6 



Km 



where the correction AE h ^ H denotes the Higgs, quark and squark mass dependent part, 
while the rest represents the total result in the limit of heavy loop particle masses. The 
latter can be derived from low-energy theorems analogously to the photonic case. The 
squark contributions to the gluonic Higgs couplings in the limit of large squark masses 
arise from the effective Lagrangian [TC = h, H] 



-I pQ( a s)/a s 
' eSf ~ *«41 + 7mo (a.) 



p V 



(23) 



where j3Q(a s )/a s = (a s /127r)[l + lla s /2n + • • •] denotes the heavy squark contribution 
to the QCD (3 function and 7 m _ (a s ) = a s /ir + ■ ■ ■ the anomalous squark mass dimension. 
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The NLO expansion of the effective Lagrangian reads a: 



AC eff = g%-^G a ^G*- 



n 



9as 

2 7T 



+ 0(a 



(24) 



Including the quark contributions in the heavy quark limit, the total effective Lagrangian 
is given by 



C 



YItx 



n 



G^-l^g* 



Q 



lla s 
1 + -r— 

4 7T . 



0* 

+ E- 

4- 4 



9 a s 
1 + — - 

2 7T . 



(25) 



The calculation based on this effective Lagrangian yields the following results for the finite 
parts of the individual contributions to the coefficients E h l H , 



E h ' H = E h J» + + E : 



-,h/H 
J 91Q 



33 - 2N F , al 
■ log '' 



6 



M 2 h/H 



E 



h/H 

virt 



71 



P]h/H 
999 



2 2 
73 



h/H A h/H ( ) 



^Q9 h Q /H A h Q /H (r Q )+EQ 9 r A r^ 



h/H ,h/H, 



E, 



h/H 
9QQ 



(26) 



which agrees with the explicit results for the coefficients E h ' H of Eq. ( J22l) in the heavy 
quark and squark limits, where AE h ^ H vanishes. 

Since at the QQ thresholds ++ states can form, the NLO QCD corrections exhibit 
Coulomb singularities as for the photonic Higgs couplings. The singular behavior can be 
derived from the Sommerfeld rescattering corrections and leads to the following expres- 
sions at each specific QoOo threshold, 



h/H A h/H 

9 Q A Qo [ " 



16tt 2 



Qo^3(7r 2 -4) 



-log(r^-i; 



ITT 



const 



^ h/H Ai/H i 
9a n n \ T C 



(27) 



which agrees quantitatively with the numerical results. 

The partial decay widths into gluons are presented in Fig. [10] with and without the 
squark contributions for two values of tg(3 = 3, 30 in the gluophobic Higgs scenario. The 
renormalization scale has been identified with the corresponding Higgs mass, /ir = M^/h- 
The partial widths range between 10" 2 and 10 MeV. The comparison of the two curves 
underlines the large size of the squark contributions, which are of the same order as the 



The value for the QCD corrections in the heavy squark limit differs from the result obtained in 
Ref. [13j . The difference can be traced back to a wrong expression for the anomalous squark mass 
dimension used in 113 . 
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quark loops in this scenario. The spikes of the full results originate from the Coulomb 
singularities of the titi,bib\ and b 2 b 2 thresholds in consecutive order. 

Fig. [11] displays the relative QCD corrections to the gluonic Higgs decays with and 
without squark contributions. Apart from the threshold singularities they are of similar 
size and increase the decay widths by 0(50%). 

In order to quantify the squark mass effects on the relative QCD corrections, Fig. [12] 
depicts the ratio between the fully massive gluonic decay widths Y{h/H — > gg) and 
the widths obtained by taking the squark contributions to the coefficients E h l H in the 
infinite squark mass limits, while the LO decay widths are used with the full squark mass 
dependence. It is clearly visible that squark mass effects on the relative NLO corrections 
can reach 20-30% in addition to the squark mass dependence at LO. Thus the fully 
massive results are significant for a proper quantitative prediction of the gluonic decay 
widths at NLO. The heavy squark mass limit turns out to be less reliable than the heavy 
top mass limit for the top quark loops. 



2.3 Gluon Fusion 



The gluon fusion processes gg — > h/H are mediated by heavy quark and squark triangle 
loops with the latter contributing significantly for squark masses below 400 GeV. The 
lowest order cross sections in the narrow-width approximation can be obtained from the 
gluonic decay widths of the scalar Higgs bosons [TfSf3"3"] . 



a LO (pp^h/H) = aQ /H T h/H 



dC" 



h/H 



71 



h/H 



dr h / H 

r LO (h/H 



(28) 



G F a 2 (fi R ) 



Q Q 



(29) 



where t^/h = M 2 j H /s with s specifying the squared hadronic cm. energy. The LO form 
factors are identical to the corresponding form factors of the gluonic decay modes, 

Eq. (J2T]) . The gluon luminosity at the factorization scale \ip is defined as 



~dV 



I — 9(x,fi 2 F )g(T/x,^ 2 F ) 

J T X 



where g(x, [i 2 F ) denotes the gluon parton density of the proton. 

The NLO QCD corrections consist of the virtual two-loop corrections, corresponding 
to the diagrams of Fig. [SI as well as the real corrections due to the radiation processes 
gg^gh/H,gq—>-qh/H and qq gh/H. The diagrams of the real corrections are shown 
in Fig. [13] While the Higgs bosons do not acquire any transverse momentum at LO, they 
appear at finite transverse momenta in these radiation processes corresponding to Higgs 
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Figure 10: QCD corrected partial decay widths of the scalar MSSM Higgs bosons to two 
gluons as functions of the corresponding Higgs masses for tg/3 = 3 and 30. The full 
curves include all loop contributions, while in the dashed lines the squark contributions are 
omitted. The kinks, bumps and spikes correspond to the titi,tt, b±bi and 6 2 &2 thresholds in 
consecutive order with rising Higgs mass. The renormalization scale of the strong coupling 
a s is chosen as the corresponding Higgs mass. 
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Figure 11: Relative QCD corrections to the partial decay widths of the scalar MSSM 
Higgs bosons to two gluons as functions of the corresponding Higgs masses for tgfl = 3 
and 30. The full curves include all loop contributions, while in the dashed lines the squark 
contributions are omitted. The kinks, bumps and spikes correspond to the t\ti,ti,bibi and 
6262 thresholds in consecutive order with rising Higgs mass. The renormalization scale of 
the strong coupling a s is chosen as the corresponding Higgs mass. 
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Figure 12: Ratio of the QCD corrected partial decay widths of the scalar MSSM Higgs 
bosons to two gluons including the full squark mass dependence and those obtained by 
taking the relative QCD corrections to the squark loops in the heavy mass limit as functions 
of the corresponding Higgs masses for tg/3 = 3 and 30. The kinks and spikes correspond 
to the titi,b\bi and b 2 t>2 thresholds in consecutive order with rising Higgs mass. The 
renormalization scale of the strong coupling a s is chosen as the corresponding Higgs mass. 
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Figure 13: Typical diagrams for the real NLO QCD corrections to the squark contributions 
to the gluon fusion processes. 



+ jet production. The relevance of quark mass effects on the shapes of the transverse 
momentum distributions has been demonstrated in Ref. [51j . so that similar effects may 
be expected for the squark mass dependence. The final results for the total hadronic cross 
sections can be split into five parts accordingly, 



a(pp -> h/H + X) = o h J H 



1 + C h/H^ 

IT 



Th/H 



dr. 



'h/H 



+ Aajf + Aajf + , (30) 



The strong coupling constant is renormalized in the MS scheme, with the top quark and 
squark contributions decoupled from the scale dependence. The quark and squark masses 
are renormalized on-shell. The parton densities are defined in the MS scheme with five 
active flavors, i.e. the top quark and the squarks are not included in the factorization 
scale dependence. The virtual coefficients split into the infrared 7r 2 term, a logarithmic 
term including the renormalization scale fiR and finite (s)quark mass dependent pieces 



C h / H (r Q , 



7T 2 + C h / H 



33 - 2N F , al 
■log 



6 



K/H 



(31) 



The finite hard contributions from gluon radiation as well as gq and qq scattering depend 
after mass factorization on the renormalization and factorization scales Hr, Hf and are 
given by 



i 

T h/H 



dC" a s h /H 
dr — - — x — o~n 



P g9 (z)l0g^ + d h g / g H (z,T Q ,TQ 



dr 



7T 
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+12 



'log(l - z) 



z[2 - z(l - z)] log(l - z) 



h/H f 1 r d£ qq CXs h/H jh/H, \ /00 n 

99 = J dT l^^^ X d qq {z,Tq,Tq), (32) 



with 2 = Th/n/T = M%, H /s, where s denotes the squared partonic cm. energy; P gg and 
P gq are the standard Altarelli-Parisi splitting functions [35] : 



P«W - ii+il^ll. (33) 



The natural scale choices turn out to be u.r = /ijr = M h / H . The quark and squark 



mass dependences are contained in the in the kernels c h ^ H (tq,Tq) and d^ H (z,tq,tq) in 

h I FT 

addition to the LO coefficients cr .In the heavy loop particle mass limit they reduce 
to simple expressions, 

h/H, \ 11 , 7 so \ Zq9q H A% H (Tq) 

c h/H (r Q ,r ) - - + -^e - . .„ — w x 




2 2 lE^r^r^+^^r^: 

d h q ^(z,T Q ,^) - ^(l-^) 3 - (34) 

These can also be derived from the effective Lagrangian Eq. (|25p . so that the full calcu- 
lation agrees with the derivation from the low-energy theorems in the heavy loop particle 
mass limits. 

Due to the formation of ++ states at the QQ thresholds the NLO QCD corrections 
develop Coulomb singularities analogous to the photonic Higgs decay modes. The singular 
behavior can be derived from the Sommerfeld rescattering corrections in the same way as 
in the photonic case, leading to the following expressions at each specific QoQ threshold, 

c m _ Se 1 9 ^ H < H{T ^^^ ,ogiT ^:! H ,r +mist ' \ , ( 35) 



e o9 wv )+e i5 4 / x / 'V< 

which agrees quantitatively with the numerical results. 
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The total gluon fusion cross section at NLO is displayed in Fig. [TH with and without 
squark loop contributions. The renormalization and factorization scales have been iden- 
tified with the corresponding Higgs mass, Hr = Hf = M h / H . In the gluophobic Higgs 
scenario the squark loops alter the size of the cross sections by up to factors of about three. 
The spikes of the full curves originate from the Coulomb singularities at the t±ti, b\bi and 
6262 thresholds in consecutive order analogously to the photonic and gluonic Higgs cou- 
plings. It can be inferred from Fig. [TH that squark effects are always important, if the 
Higgs mass exceeds the corresponding squark- ant isquark threshold. This is confirmed in 
particular for large values of tg/3, where the top quark contribution is less important, 
but sizeable squark effects are visible close to and beyond the different thresholds. This 
feature also holds in other MSSM scenarios. 

The LO and NLO cross sections are shown in Fig. [15j The QCD corrections increase 
the gluon fusion cross sections by 10-100%, but can be significantly larger in regions 
of large destructive interferences between quark and squark loops, as is the case for very 
large Higgs masses for tg/3 = 30. The corrections are of very similar size for the quark and 
squark loops individually in agreement with the results of Ref. [13]. In spite of the large 
corrections the residual scale dependence is reduced from about 50% at LO to ~ 20% at 
NLO and indicates a significant stabilization of the theoretical predictions. This agrees 
with the former results for the quark loop contributions [8fT0] . Based on the approximate 
NNLO and NNNLO results in the limit of heavy top quarks a further moderate increase 
by less than 20-30% can be expected beyond NLO for small and moderate values of tg/3. 
For large values of tg/3, however, the size of the NLO corrections is moderate in regions 
without strong destructive interference effects between the quark and squark loops and 
the scale dependence is small. This signalizes a much more reliable result after including 
the NLO corrections. The residual theoretical uncertainties of our NLO results can be 
estimated to less than about 20%. 

The squark mass effects on the K factors are exemplified in Fig. [TBI where the ratios 
of the NLO cross sections are displayed, including the full mass dependence, and of the 
NLO cross sections, where the coefficients c h ^ h (tq, ta) and d^j H (z, tq, tq) are used in the 
heavy squark limits (tq — > 00). In addition to the squark mass dependence of the LO 
cross sections, the K factors develop a squark mass dependence of up to about 20%, thus 
supporting the relevance of our results compared to the previous results of Ref. [13] . The 
squark mass effects on the K factors turn out to be larger than the corresponding quark 
mass effects [9]. In addition they are larger than the residual theoretical uncertainties and 
cannot be neglected in realistic analyses. Since the gluino contributions are expected to 
be much smaller, the squark mass dependence obtained in this work will be the dominant 
part of the differences between the heavy mass limits and a full MSSM calculation at 
NLO. 
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Figure 14: QCD corrected production cross sections of the scalar MSSM Higgs bosons via 
gluon fusion as functions of the corresponding Higgs masses for tg/3 = 3 and 30. The full 
curves include all loop contributions, while in the dashed lines the squark contributions are 
omitted. The kinks, bumps and spikes correspond to the titi,ti,bibi and b<^>2 thresholds 
in consecutive order with rising Higgs mass. The renormalization and factorization scales 
are chosen as the corresponding Higgs mass. 
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Figure 15: Production cross sections of the scalar MSSM Higgs bosons via gluon fusion as 
functions of the corresponding Higgs masses for tg/3 = 3 and 30. The full curves include 
the QCD corrections, while the dashed lines correspond to the leading-order predictions. 
The kinks and spikes correspond to the and thresholds in consecutive order 

with rising Higgs mass. The renormalization and factorization scales are chosen as the 
corresponding Higgs mass. 
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Figure 16: Ratio of the QCD corrected production cross sections of the scalar MSSM Higgs 
bosons via gluon fusion including the full squark mass dependence and those obtained by 
taking the relative QCD corrections to the squark loops in the heavy mass limit as functions 
of the corresponding Higgs masses for tg/3 = 3 and 30. The kinks and spikes correspond 
to the titi,bibi and 6 2 ^2 thresholds in consecutive order with rising Higgs mass. The 
renormalization and factorization scales are chosen as the corresponding Higgs mass. 
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3 Conclusions 



We have presented a NLO QCD calculation for the squark loop contributions to the pro- 
duction of neutral scalar MSSM Higgs bosons at the LHC and their decay modes into 
gluons and photons. The photonic couplings of these Higgs particles play a significant role 
at a future photon collider, which may be built by Compton backscattering of laser light 
from electron beams of a linear e + e~ collider. The corrections stabilize the theoretical 
predictions compared to the LO predictions. The QCD corrections turn out to be size- 
able for the photonic Higgs couplings and large for the gluonic Higgs decays as well as the 
gluon fusion processes. They increase the latter cross sections and gluonic decay widths 
significantly so that the results have to be taken into account for reliable analyses based on 
these processes. Despite of the large size of the NLO corrections the approximate results 
in the heavy top mass limits for the gluonic decay widths and gluon-fusion cross sections 
indicate sufficient perturbative convergence. The corrections beyond NLO are expected 
to be of moderate size for the central scale choices. Squark mass effects on the relative 
QCD corrections are sizeable and larger than the corresponding quark mass effects for 
the quark loop contributions. It should be noted that squark mass effects are always 
relevant, if the squark masses are of the order of the top mass, or the Higgs mass is larger 
than the corresponding virtual squark-antisquark threshold. Our results have been imple- 
mented in the program HIGLU [36J and can thus be applied to other MSSM scenarios, too. 

Note added in proof. During the final write-up of our work two independent papers 
appeared [37,38J, where the virtual corrections to the quark and squark loops of the gluon 
fusion processes gg — > h/H have been derived analytically. However, a full numerical 
analysis of the gluon fusion processes at NLO is not contained in these papers. We have 
compared our virtual corrections with the authors of Ref. [38] and found full numerical 
agreement for the virtual corrections to the photonic and gluonic Higgs couplings. 
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